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Bond stretching mimics different levels of electron correlation and provides a challenging testbed 
for approximate many-body computational methods. Using the recently developed phaseless 
auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well- 
studied molecules BH and Na, and in the H50 chain. To control the sign/phase problem, the 
phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an ap- 
proximate phase condition that depends on a trial wave function. With single Slater determinants 
from unrestricted Hartree-Fock (UHF) as trial wave function, the phaseless AF QMC method gen- 
erally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) 
method in mapping the potential-energy curve. In both BH and N2, we also study the use of 
multiple-determinant trial wave functions from multi-configuration self-consistent-field (MCSCF) 
calculations. The increase in computational cost versus the gain in statistical and systematic accu- 
racy are examined. With such trial wave functions, excellent results are obtained across the entire 
region between equilibrium and the dissociation limit. 
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I. INTRODUCTION 

Quantum Monte Carlo (QMC) methods are an at- 
tractive means to treat explicitly the interacting many 
fcrmion system. Their computational cost scales favor- 
ably with system size, as a low power. The ground 
state wave function is obtained stochastically by Monte 
Carlo (MC) sampling, either in particle coordinate space 
P, or in Slater determinant space [1, 0|. Except for 
a few special cases, however, these methods suffer from 
the fermion sign problem 0, [f|, which, if uncontrolled, 
causes an exponential loss of the MC signal and negates 
the favorable computational scaling. No formal solution 
has been found for this problem, but approximate meth- 
ods have been developed that control it. These include 
the fixed-node method Q in real coordinate space and 
constrained path methods 0, 0, @| in Slater determi- 
nant space. The real-space fixed-node diffusion Monte 
Carlo (DMC) method has long been applied to a vari- 
ety of solids and molecules [l[. Recently, the phaseless 
auxiliary-field (AF) method was introduced which pro- 
vides a framework for ab initio electronic structure cal- 
culations by QMC in Slater determinant space, within a 
Hilbert space defined by any single-particle basis 0]. 

The phaseless AF QMC method controls the phase 
problem in an approximate way by using an input trial 
wave function (WF) [|. This is a generalization of the 
constrained path approach Q which has been applied to 
lattice models with Hubbard-like interactions. Compared 
with previous efforts [1, Q on realistic electronic systems 
using the standard auxiliary-field formalism [Tol . flTI | , the 
phaseless AF QMC method overcomes the poor (expo- 
nential) scaling with system size and projection time and 
has statistical errors that are well-behaved. 
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The systematic error from the phaseless approximation 
has been found to be small near equilibrium geometries 
in a variety of systems. The method was applied using 
a planewave basis with pseudopotentials to several sp- 
bonded atoms, molecules, and solids [!, [H, [lj| and to 
the transition metal molecules TiO and MnO [l4j]. It 
has also been applied, with Gaussian basis sets, to first- 
and second-row atoms and molecules [Hj], to post-d el- 
ements (Ga-Br and In-I) and to hydrogen bonded 
systems [ljj]- The calculated all-electron total energies 
of first-row atoms and molecules at equilibrium geome- 
tries show typical systematic errors of no more than a 
few milli-hartree (mE n ) compared to exact results. This 
accuracy is roughly comparable to that of CCSD(T), cou- 
pled cluster with single- and double-excited clusters plus 
a non-iterative correction to the energy due to triple ex- 
cited clusters. In post-d systems, our results with several 
basis sets are in good agreement with CCSD(T) results 
and, for l arge basis sets, in excellent agreement with ex- 
periment In almost all of these calculations, we have 
used as trial WF mean-field solutions from independent- 
electron calculations. 

Bond stretching provides a difficult test for approxi- 
mate correlated methods. In the dissociation limit, the 
unrestricted Hartree-Fock (UHF) solution gives a qual- 
itatively correct description of the system. The inter- 
mediate region between the equilibrium and dissociated 
geometries represents a situation analogous to a metal- 
insulator transition. Due to quasi-degeneracies, there 
can be more than one important electronic configura- 
tion, and a single determinant often cannot adequately 
describe the system. Multi-configurational approaches 
can describe to a large degree the static correlations in 
the system, but often miss a large proportion of the dy- 
namic correlations. 

No general method has demonstrated the ability to 
consistently maintain uniformly high accuracy away from 
equilibrium. Coupled cluster (CC) methods [I&llSj], sucu 
as CCSD(T), are remarkably good in describing the equi- 
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librium properties, but are less successful in describing 
systems with quasi-degeneracies such as the case in the 
breaking of chemical bonds (20I . I21I I22I I23I ] . Higher or- 
der clusters have to be fully included in the iterative ap- 
proach, because the perturbative corrections are based on 
non-degenerate perturbation theories, and usually lead 
to divergences for stretched nuclear geometries. Since 
CCSD already scales as TV 6 with basis size, going to triple 
and higher order clusters is computationally expensive. 
Multi-reference CC methods could potentially solve some 
of these problems, but unlike the single-reference CC 
method, these are still not widely established [H| . Other 
coupled-cluster-based approaches have been introduced 
recently to handle bond stretching, and this remains an 
active field of research; see for example Refs. [25L [26l. [27lj . 

In this paper, we test the phaseless AF QMC method 
away from Born-Oppenheimer equilibrium configura- 
tions. We investigate bond stretching in two well- 
studied molecules, BH and N 2 , and in a hydrogen chain, 
H50, where exact or very accurate results from full- 
configuration interaction ( FCI) or density-matrix renor- 
malization group (DMRG) [H, [HHfil are available. We 
first use single Slater determinant trial WF's, obtained 
by the unrestricted Hartree-Fock (UHF) method. It is 
shown that AF QMC with UHF as trial WF generally 
gives better overall accuracy and a more uniform be- 
havior than CCSD(T). The use of multiple determinant 
trial WFs from multi-configuration self-consistent-field 
(MCSCF) calculations is then examined in the diatomic 
molecules. With these trial WFs, excellent AF QMC 
results are achieved across the entire potential energy 
surface. 

The rest of the paper is organized as follows. The 
phaseless AF QMC method is first briefly reviewed in 
the next section. In Sec. IIII1 we present and discuss the 
potential-energy curves of the various systems. Finally, 
in Sec. IIV[ we conclude with a brief summary. 

II. THE PHASELESS AF QMC METHOD 

The many-body Born-Oppcnhcimcr Hamiltonian in 
electronic systems can be written in second quantization, 
in any single-particle basis, as 

N N 

H = H 1+ H 2 = J2 T ^i°i + 2 E V mi44c kCl , (1) 

i,3 i,j,k,l 

where N is the size of the chosen one-particle basis, and 
c\ and Ci are the corresponding creation and annihilation 
operators. The one-electron Xy and two-electron Vijki 
matrix elements depend on the chosen basis. 

The phaseless AF QMC obtains the ground state of 
the system by projecting from a trial WF \^t) which 
has a non zero overlap with the exact ground state of the 
system: 

/ - \ M 

\* GS ) =Vm o (e- TH ) |* r ), (2) 



where r is a small time-step, and \^t) is assumed 
to be in the form of a Slater determinant or a lin- 
ear combination of Slater determinants. Using a sec- 
ond order Trotter decomposition, we can write e~ rH = 
e -rH 1 /2 e -rH 2e -TH 1 /2^ Tne resmtmg Trotter time-step 
error decreases with t, and can be eliminated by an ex- 
trapolation to t — with multiple calculations. 

The central idea in the AF QMC method is the use of 
the Hubbard-Stratonovich (HS) transformation [3l| : 

e-^=n(4= [°° dv a e-rte^°-V£*-). (3) 

a \v27Tj_oo J 

Equation ([3]) introduces one-body operators v a which can 
be defined generally for any two-body operator by writing 
the latter in a quadratic form, such as H2 = — \ Y] a CaV^, 
with £ Q a real number. The many-body problem as de- 
fined by H2 is now mapped onto a linear combination 
of non-interacting problems defined by v a , interacting 
with external auxiliary fields. Averaging over different 
auxiliary-field configurations is then performed by MC 
techniques. Formally, this leads to a representation of 
|^Gs) as a linear combination of an ensemble of Slater 
determinants, { \4>) }. The orbitals of each \4>) are written 
in terms of the chosen one-particle basis and stochasti- 
cally evolve in imaginary time. 

Generally, the AF QMC method suffers from the sign 
or phase problem 0, [|| . The phaseless AF QMC method 
[H used in this paper controls the phase/sign prob- 
lem in an approximate manner using a trial WF, \^t)- 
The method recasts the imaginary-time path integral as 
a branching random walk in Slater-determinant space 
d, HI- It uses the overlap {^ T \cj)), to construct phase- 
less random walkers, \(j>) / (^xl'P) : which are invariant un- 
der a phase gauge transformation. The resulting two- 
dimensional diffusion process in the complex plane of the 
overlap ("J/tI^) is then approximated as a diffusion pro- 
cess in one dimension. The ground-state energy com- 
puted with the so-called mixed estimate is approximate 
and not variational in the phaseless method. The error 
depends on the quality of | 'J't) , and the method becomes 
exact as the trial WF approaches the exact ground state 
of the system. This is the only error in the method that 
cannot be eliminated systematically. 

In most applications to date @, [ll, [H 0, [H [H, El , 
the trial WF has been a single Slater determinant taken 
directly from mean-field calculations. We have found [HI, 
that using the UHF solution leads to better QMC 
energies than using the restricted Hartree-Fock (RHF) 
Slater determinant. This was the case even with singlets. 

In this study, we will present, in addition to the 
single-determinant trial WF, results based on multi- 
determinant trial WFs obtained from MCSCF calcula- 
tions. In some cases, such as bond stretching, a multi- 
determinant trial WF can capture some of the static cor- 
relation in the system, and thus improve the quality of 
the constraint in the phaseless approximation. A better 
trial WF will generally reduce the systematic errors of 
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FIG. 1: Potential-energy curves of the BH molecule from 
FCI, coupled cluster, and QMC methods, using the cc-pVDZ 
basis. The QMC/UHF and QMC/MCSCF curves are ob- 
tained respectively with single determinant UHF and multi- 
determinant truncated MCSCF trial wave functions. The in- 
set shows the deviations (in mEh) of the various methods from 
the FCI results. 




# determinants 

FIG. 2: Phaseless AF QMC systematic errors versus the num- 
ber of determinants included in the trial WF from MCSCF. 
The discrepancy between QMC and FCI energies in BH/cc- 
pVDZ is shown for geometries of 2R e , 2.5 R e , and 5-R e (in 
mEh). The corresponding errors from the single-determinant 
UHF trial WF are also shown, as symbols on the left. 



the phaseless AF QMC method. 

In addition, a better trial WF will typically also lead 
to better statistics in the AF QMC method, for a fixed 
number of independent MC samples. A simple measure 
of the efficiency of the multi-determinant MCSCF trial 
WF relative to a single-determinant UHF trial WF is the 
following ratio, rj: 

_ (^sample £ 2 ) MCSCF 

(AT ' ^ ' 

(/'sample e Jrjjjp 

where e is the final statistical error, and N sa _ mp i c is the 
total number of MC samples used in the calculation. (A 
more precise but closely related measure is the ratio of 
the variances of the local energy. For its purpose here as 
a rough indicator, however, the difference between them 
is not significant.) We expect r] < 1 for a reasonable 
number of determinants in the MCSCF; in general, the 
better \^t), the smaller 77. Since the computational cost 
of the phaseless AF QMC method increases linearly with 
the number of determinants in I^t), the overall com- 
putational cost of the QMC calculation with respect to 
the single-determinant trial WF is 77 times the number of 
determinants in j^r)- 

III. RESULTS AND DISCUSSION 

To examine the performance of the phaseless AF 
QMC method in bond stretching, the potential-energy 
curves of the diatomic molecules BH and N2 are 
first studied and compared to exact FCI, near-exact 
DMRG, and several levels of CC methods. In addi- 
tion, the symmetric and asymmetric bond stretching 
of an H50 linear chain is examined and compared to 
DMRG results. For the diatomic molecules, both single- 
determinant UHF (QMC/UHF) and multi-determinant 



MCSCF (QMC/MCSCF) trial WFs were used. AF QMC 
for the H50 chain used single-determinant Hartree-Fock 
trial wave functions. 

In our calculations, all the electrons are correlated and 
the spherical harmonic (as opposed to Cartesian) form of 
the Gaussian basis functions was used. For the molecules, 
cc-pVDZ basis sets were used, except in the challenging 
case of the (triple) bond stretching of N 2 , where calcu- 
lations were also performed with the cc-pVTZ basis set 
[32| . For the H50 chain, the minimal STO-6G basis set 
was used. 

All of the Hartree-Fock, MCSCF, and CC calculations 
were carried out using NWCHEM [33[ within C21, sym- 
metry. Some of these calculations were also verified using 
Gaussian 98 H and MOLPRO [35]. The MCSCF en- 
ergies were obtained from a complete-active-space SCF 
(CASSCF) [36| calculations. In most of the molecules, we 
used the RHF and UHF reference states for the CC calcu- 
lations [e.g., labeled RCCSD(T) and UCCSD(T), respec- 
tively]. FCI calculations were performed using MOLPRO 

A. BH 

Table U summarizes the cc-pVDZ basis set energies [in 
hartrees (Eh)] obtained with a variety of methods for 
seven BH geometries over a range R/R e = 1 — 5, where 
R e = 1.2344 A. The MCSCF energy was obtained by a 
CASSCF calculation, performed with 4 active electrons 
and 8 active orbitals. FigureQ]shows the potential-energy 
curves from selected methods. 

Near the equilibrium geometry, the RCCSD(T) ener- 
gies are in good agreement with the FCI energy. How- 
ever, this agreement deteriorates for larger nuclear sep- 
aration, and RCCSD(T) shows an unphysical dip for 
R > 2.5i? e which increases for larger bondlength R. The 
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TABLE I: All-electron total energies of BH versus bondlength as calculated by a variety of methods, using a cc-pVDZ basis set. 
The exact results are given by FCI. AF QMC energies obtained with both the UHF and MCSCF (see text) trial wave functions 
are shown. All energies are in Eh and statistical errors in QMC (shown in parentheses) are on the last digit. An equilibrium 
bond length of R e = 1.2344 A is used. 
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1 


1.5 


2 


2.5 


3 


4 


5 


RHF 






—94 qqo 7c;q 


— 94 q4fl 9^fi 


— 94 qfl9 889 


— 94 S^q 1 1 n 


— 94 84fl 0P,<\ 


UHF 


-25.131559 


-25.065 817 


-25.034 695 


-25.030 226 


-25.029 455 


-25.029 262 


-25.029 241 


MCSCF 


-25.199 413 


-25.150 182 


-25.105 151 


-25.086 687 


-25.081836 


-25.080 626 


-25.080 571 


RCCSD 


-25.214 360 


-25.163 701 


-25.112 404 


-25.088 577 


-25.080 907 


-25.078 318 


-25.078 097 


RCCSD(T) 


-25.215 767 


-25.165 880 


-25.117 034 


-25.099 308 


-25.100152 


-25.107 421 


-25.109 133 


UCCSD 


-25.214 360 


-25.163 498 


-25.110 239 


-25.091 729 


-25.089 185 


-25.088 720 


-25.088 710 


UCCSD(T) 


-25.215 767 


-25.165 784 


-25.114 131 


-25.093 105 


-25.090 055 


-25.089 555 


-25.089 545 


QMC/UHF 


-25.214 9(3) 


-25.164 0(2) 


-25.114 0(1) 


-25.093 9(3) 


-25.088 7(1) 


-25.086 9(4) 


-25.0871(2) 


QMC/MCSCF 


-25.215 89(6) 


-25.166 51(7) 


-25.117 78(8) 


-25.097 97(6) 


-25.092 74(9) 


-25.091 58(6) 


-25.091 3(2) 


FCI 


-25.216 249 


-25.166 561 


-25.117 705 


-25.097 084 


-25.091467 


-25.089 986 


-25.089 912 



failure of RCCSD (T) to describe the molecule for larger 
bondlengths is attributed to the poor quality of the RHF 
WF in describing bond breaking. In the large bondlength 
limit, the UHF solution is better than the RHF solution. 
This is reflected also in CC results based on the UHF 
solution; the UCCSD(T) energies are in very good agree- 
ment with the FCI energy for large R. The UCCSD(T) 
energies are in less good agreement with FCI in the inter- 
mediate region. Overall, UCCSD (T) does quite well in 
BH, which has a relatively small number of excitations. 

As shown in Table |H QMC/UHF energies are com- 
parable to RCCSD(T) and in good agreement with FCI 
near the equilibrium geometry. As the bond is stretched, 
QMC/UHF energies become less accurate. The discrep- 
ancy with FCI energies is ~ 3 rnEi, for R > 2 R e . In the 
QMC/MCSCF calculations, the multi-determinant trial 
WF included determinants from MCSCF with coefficient 
cut-offs > 0.01. Thus the variational energy of our MC- 
SCF WF is higher than the corresponding MCSCF result 
listed in the tables. The average value of rj, as defined in 
Eq. (U]), is 0.04, and the largest value is 0.08, at the largest 
bondlength. The QMC/MCSCF energies are in excellent 
agreement with the FCI energies, to within w 1 mEh for 
all studied bondlengths. 

The optimum cutoff value of determinant coefficient 
cut-off in the MCSCF trial WF is, of course, system- 
dependent. The accuracy of the QMC calculation gener- 
ally improves as the cutoff is lowered, while the compu- 
tational cost increases. For a small system like BH, a rel- 
atively low cutoff leads to excellent trial WFs with large 
efficiency gain, as the r\ values show. Figure [5] shows the 
QMC errors as a function of the number of determinants 
included in the trial WF for three geometries of BH. For 
2 R e , the QMC results with MCSCF trial WFs containing 
determinants with coefficient cut-offs less than 0.02 (29 
determinants) and less than 0.01 (52 determinants) are 
equivalent within statistical errors. Similarly, for 2.5 i? e , 
the QMC results obtained with trial WFs of 24 and 44 
determinants are indistinguishable within the statistical 
errors. Indeed, in both of these cases, 8 determinants 
in the trial WF give systematic errors less than 2 mEh. 



By contrast, for 5i? e , considerably more determinants 
are required to achieve converged QMC systematic er- 
rors. Note that, because the MCSCF WF is in a spin 
restricted form, more than one determinant (many more 
in the case of 5 R e ) is required to surpass the accuracy of 
QMC/UHF. 

B. N 2 

Bond stretching in N2 is particularly challenging, be- 
cause it involves the breaking of a triple bond. As a re- 
sult, N2 has been extensively studied |22l. I23I l38l. l39l |4C| . 
Table |TT] summarizes the calculated total energies, us- 
ing cc-pVDZ and cc-pVTZ basis sets. Figure [3] plots a 
selected subset of these potential-energy curves. With 
the cc-pVDZ basis set, CC results based on the UHF 
reference state, and the near-exact DMRG energies are 
from Ref. @. We have also verified the UCCSD and 
UCCSD(T) energies. For both basis sets, the CASSCF 
calculations are performed with 6 active electrons and 12 
active orbitals. 

The main features of the CC potential-energy curves 
of N2 are similar to those of BH. In contrast with the BH 
molecule, however, the effects beyond double excitations 
are substantial in N2, even at the equilibrium geometry. 
CC results based on a RHF reference show an unphysi- 
cal dip for R > 3.6 bohr (R/R e > 1.75 in Fig. ©. For 
the cc-pVDZ basis at the larger R = 3.6 and 4.2 bohr 
bond lengths, UCCSD (T) based on the UHF reference is 
in better agreement with DMRG than RCCSD(T). Fully 
including triple excitations with UCCSDT leads to a sig- 
nificant improvement over UCCSD(T) for all geometries 
except R e , while RCCSDT seems to be slightly worse 
than RCCSD (T), except at the last geometry. 

QMC with an UHF trial wave function gives a bet- 
ter overall accuracy and a more uniform behavior than 
CCSD(T) in mapping the potential-energy curve in the 
cc-pVDZ basis. The largest difference of the QMC/UHF 
energies compared to DMRG is at the second to last nu- 
clear separation, and is approximately 9 mEh- With 
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TABLE II: All-electron energies of N2 versus bondlength (in bohr) as calculated by a variety of methods, using cc-pVDZ 
and cc-pVTZ basis sets. With the cc-pVDZ basis, DMRG energies and CC results based on a UHF reference are taken from 
Ref. QMC/UHF energies are obtained with the UHF trial WF, and QMC/MCSCF are obtained with a multi-determinant 
trial WF from an MCSCF calculation (see text). All energies are in Eh, and QMC statistical errors (shown in parentheses) are 
on the last digit. 
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QMC/MCSCF, we included in the multi-determinant 
trial WF all determinants with a weight larger than 0.01. 
This gives 65, 66, 76, 97, 82, and 58 determinants for 
the six bondlengths (in ascending order), respectively. 
As can be seen from Table [TT] and the inset of Fig. [31 
the agreement between the QMC/MCSCF and DMRG 
values is more uniform and the discrepancy is less than 
2-3 mEh for all geometries. 

In the QMC/MCSCF calculations for the cc-pVDZ ba- 
sis set, the average value of rj of Eq. ^ is 0.42, and the 
largest value is 0.80 at the equilibrium geometry. The 
weight cutoff choice of 0.01 in selecting the determinants 
to include from the MCSCF WF was the same as in the 
BH calculations. This was likely too conservative as in 
BH. For example, with R = 2.118 bohr, the QMC results 
were within statistical errors for a trial WF that included 
determinants with coefficient cut-offs > 0.035. 

The cc-pVTZ results from the various methods parallel 
very well the cc-pVDZ results, as can be seen from Fig. [3] 
and Table UU Both the QMC/UHF and QMC/MCSCF, 
for example, mirror each other in the two basis sets. We 
thus expect the accuracy of the different QMC and CC 
methods using the cc-pVTZ basis to be comparable to 
that using the cc-pVDZ basis, where DMRG results are 
available. For the cc-pVTZ basis set, the QMC/MCSCF 



calculations included determinants with coefficient cut- 
offs > 0.02. The average value of 77 is 0.16 and the 
largest value is 0.41 for R = 2.7 bohr. Additional 
QMC/MCSCF calculations were performed for R = 4.2, 
3.6, and 2.7 bohr, including determinants with coefficient 
cut-offs > 0.01, and the same energies were obtained as 
those in Table HT1 within statistical errors. 



C. Hydrogen chain: H50 

The hydrogen linear chain exhibits characteristic sig- 
natures of a metal-insulator transition as the interatomic 
distances are varied. It also provides a simple but chal- 
lenging model for extended systems, where the favorable 
scaling of QMC will be especially valuable. Bond stretch- 
ing in a linear chain of hydrogen atoms, H50, was recently 
benchmarked with DMRG [41( . This 50-electron system 
was treated using a minimal STO-6G basis set of 50 or- 
bitals. Both symmetric and asymmetric bond stretching 
were considered. In the case of symmetric bond stretch- 
ing, the bond between consecutive hydrogen atoms is 
stretched over the range R = 1.0 — 4.2 bohr, and the final 
structure consists of 50 equidistant, nearly- independent 
H-atoms. In the case of asymmetric bond stretching, 25 
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TABLE III: Symmetric and asymmetric bond stretching in a H50 linear chain, using a minimal STO-6G basis set. Total energies 
versus geometry are shown for different methods. DMRG and RCCSD(T) values are from Ref. [4l|]. Bondlengths are in bohr 
and energies are in Eh. QMC statistical errors (shown in parentheses) are on the last digit. 



R 


RHF 


UHF 


RCCSD(T) 


DMRG 


AFQMC 


Symmetric 












1.0 


-16.864 88 


-16.864 88 


-17.282 27 


-17.284 07 


-17.285 2(1) 
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-22.46127 


-22.468 05 


-22.944 57 


-22.947 65 


-22.947 5(7) 
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— ZO.Uzy / 


— zo.uoo yi 
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1.6 


— 26.062 25 


-26.130 19 


-26.713 14 


-26.719 44 


-26.718 8(5) 


1.8 


— 26.265 98 


-26.396 69 


-27.03145 


-27.038 65 


-27.038 8(3) 


2.0 


-26.008 20 


-26.237 77 


-26.920 90 


-26.926 09 


-26.925 6(9) 


2.4 


-24.835 76 


-25.434 02 




-26.160 57 


-26.159 4(5) 


2.8 


-23.360 81 


-24.63419 




-25.274 80 


-25.276 5(7) 


3.2 


-21.896 33 


-24.108 60 




-24.568 28 


-24.573 3(5) 


3.6 


-20.574 29 


-23.823 26 




-24.102 77 


-24.108 4(7) 


4.2 


-18.955 95 


-23.634 41 




-23.749 71 


-23.748 9(4) 


-rvoy iiiiiit; LI 1L, 












1.6 


-25.963 71 




-26.486 01 


-26.487 38 


-26.486 7(7) 


1.8 


-26.61768 




-27.126 41 


-27.12716 


-27.126 4(9) 


2.0 


-27.07182 




-27.576 79 


-27.577 32 


-27.576 9(3) 


2.4 


-27.609 24 




-28.11727 


-28.117 61 


-28.116 3(7) 


2.8 


-27.873 62 




-28.386 84 


-28.387 07 


-28.384 2(5) 


3.2 


-28.004 68 




-28.521 10 


-28.52124 


-28.518 9(7) 


3.6 


-28.069 65 




-28.587 28 


-28.58736 


-28.584 4(5) 


4.2 


-28.11100 




-28.628 54 


-28.628 58 


-28.626 5(5) 



equivalent H 2 molecules are considered, each with a fixed 
bondlength of 1.4 bohr, where two consecutive hydrogen 
atoms belonging to two different H 2 molecules are sepa- 
rated over a range of R — 1.4 — 4.2 bohr, and the final 
structure consists of 25 equidistant, nearly- independent 
H2 molecules, each at its equilibrium bondlength. Ta- 
ble IIIII shows the results for both symmetric and asym- 
metric bond stretching. The RHF and UHF energies, as 
well as our QMC results obtained using the UHF trial 
wave function (or RHF when there is no UHF solution) 
are shown. The RCCSD(T) and DMRG energies as re- 
ported in Ref. [4l| are also shown for comparison. 

For symmetric stretching, the QMC/UHF total ener- 
gies are in good agreement with the DMRG results, with 
the largest discrepancy being about 5 mEh for R = 3.2 
and R = 3.6 bohr. As the bondlength is stretched, 
the correlation energy of the system increases. In view 
of the above results for bond stretching in diatomic 
molecules, it is not surprising that the discrepancy be- 
tween RCCSD(T) and DMRG increases as R is increased, 
and for R > 2.4, RCCSD(T) fails to converge as reported 
in Ref. 

For asymmetric bond stretching, the QMC energies are 
again in good agreement with the DMRG values. The 
difference between the QMC and DMRG total energies 
is less than as 2 — 3 mEh for all bond lengths. Here 
no distinct UHF solution was found, so the RHF Slater 
determinant was used as the trial WF. The RHF trial 
WF dissociates properly as R is increased in this case, so, 
not surprisingly, RCCSD(T) is in good agreement with 
DMRG. 



IV. SUMMARY 



Bond stretching in chemistry is a non-trivial challenge 
for all approximate correlated methods. In this paper, we 
applied the recently introduced phaseless auxiliary-field 
QMC method to study bond-stretching in BH and N2 
and in the H50 chain. The quality of the phaseless AF 
QMC method depends on the trial wave function that 
is used to control the sign/phase problem. With a single 
UHF Slater determinant as trial wave function, AF QMC 
has performed very well for molecular geometries near the 
equilibrium configuration, as shown by comparisons with 
exact values, CCSD(T) calculations, and experimental 
results. The results in this paper are consistent with this 
and extend the reach of phaseless AF QMC method be- 
yond Born Oppcnheimer equilibrium structures to bond 
stretching and bond breaking. For larger nuclear separa- 
tions, we find that AF QMC with a single-determinant 
UHF solution in general gives better overall accuracy and 
a more uniform behavior than coupled cluster CCSD(T). 
In some stretched bond situations, however, QMC/UHF 
errors are seen to be significant. In these cases, we find 
that a trial wave function with a modest number of de- 
terminants usually reduces the QMC error to a few mEh ■ 
The QMC computational cost of the multi-determinant 
trial wave function scales linearly with the number of 
determinants, but a better trial wave function can re- 
duce both systematic and statistical errors. Using multi- 
determinant trial wave functions taken directly from MC- 
SCF calculations, the AF QMC results are in very good 
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FIG. 3: Potential-energy curves of Na, using cc-pVDZ and 
cc-pVTZ basis sets. QMC/UHF energies are obtained with a 
UHF trial WF, and QMC/MCSCF with a truncated multi- 
determinant trial WF taken from a MCSCF calculation. CC 
results at several levels are shown for both RHF and UHF 
reference states. For the cc-pVDZ basis set, DMRG results 
from Ref. [3£| are also shown, and the inset shows the de- 
viations (in mEh) of the various methods from the DMRG 
potential-energy curve. 
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